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Abstract. Recent modeling of multi- waveband spectroscopic and maser observations suggests 
that the ionized outflows in the nuclear region of the archetypal Seyfert-2 galaxy NGC 1068 
are inclined with respect to the vertical axis of the obscuring torus. Based on this suggestion, 
we build a complex reprocessing model of NGC 1068 for the optical/UV band. We apply the 
radiative transfer code STOKES to compute polarization spectra and images. The effects of 
electron and dust scattering and the radiative coupling occurring in the inner regions of the 
multi-component object are taken into account and evaluated at different polar and azimuthal 
viewing angles. The observed type-l/type-2 polarization dichotomy of active galactic nuclei 
is reproduced. At the assumed observer's inclination toward NGC 1068, the polarization is 
dominated by scattering in the polar outflows and therefore it indicates their tilting angle with 
respect to the torus axis. While a detailed analysis of our model results is still in progress, we 
briefly discuss how they relate to existing polarization observations of NGC 1068. 



1. Introduction 

According to the unified model of active galactic nuclei (AGN) [1] the accreting supermassive 
black hole of so-called type-2 thermal AGN lies hidden behind an optically thick, dusty torus. 
The obscuration prevents the broad line region to be seen and therefore type-2 AGN lack 
the kineticly broadened Balmer emission lines that characterize type-1 objects. Observations 
indicate that outflows are ejected from close to the accretion disk and it is possible that 
the torus funnel collimates these winds along polar directions. The launch of the wind is 
not yet understood, but it is a natural assumption that the ejection and collimation happen 
symmetrically with respect to the axis of the dusty torus and the accretion disk. In the case of 
the archetypical Seyfert-2 galaxy NGC 1068, the AGN is observed at a rather high inclination 
and the ionized winds are therefore well seen from the side. 

A recent study by [2j suggests that the polar winds of NGC 1068 are inclined with respect to 
the torus axis. Based on this assumption, [3] and [J| started to investigate the radiative transfer 
effects of the misaligned outflows on the induced X-ray polarization. The modeling suggests 
that the position angle of the wind on the sky plane could be determined unambiguously from 
a soft X-ray polarization measurement. In addition to this, the tilting angle of the wind with 
respect to the torus axis could be constrained by the rotation of the polarization position angle 
between the soft and hard X-ray bands. 



flared disk 


dusty torus 


polar outflows 


R min = 0.02 pc 

flmax = 0.04 pc 

half-opening angle = 20° 
equat. optical depth = 1 
electron scattering 


-Rmin = 0.1 PC 
-Rmax = 0.5 PC 

half-opening angle = 60° 
equat. optical depth = 750 
Mie scattering 


-Rmin = 0.3 PC 
-Rmax = 1.8 PC 

half-opening angle = 40°, tilted by 18° 
vertical optical depth = 0.03 
electron scattering 



Table 1. Parameters of the three scattering regions. The polar outflows and torus half-opening 
angles are measured with respect to the vertical symmetry axis of the torus. The half-opening 
angle of the flared disk is measured from the equatorial plane. 



The predictions for future X-ray polarization measurements are more solid when they can be 
coherently related to present-day polarization results from a different waveband. Therefore, we 
here present some results for the optical/UV counterpart of the same model setup. Similarly to 
[3J, we present spectral and polarization modeling of NGC 1068 at several polar and azimuthal 
viewing angles. We use the new polarization mapping capacity that was recently implemented 
in our code stokes. A more detailed description of the new stokes version and more results 
of the current study will be published elsewhere (Marin et al., in preparation). 

2. The model setup for NGC 1068 

This research note focuses on the expected polarization signatures of NGC 1068 with misaligned 
outflows as a function of the polar and azimuthal viewing angle. We base our model on the 
description given in [2 J assuming that the bi-conical outflows are inclined by 18° from the vertical 
torus axis. The model setup is very similar to the one used in Its main characteristics are 
summarized in Table [U Note that here we did not include an irradiated accretion disk as it was 
done for the X-ray band where a lamp-post geometry was assumed. Instead, the negligibly small 
accretion disk produces an optical/UV continuum spectrum having a F*(E) oc E~ a power law 
shape. The spectral index is set to a = 1. The lower and upper limits of the spectrum are set to 
2000 A and 8000 A, respectively. The Cartesian coordinate system of the model space and its 
relation to the polar and azimuthal viewing angle is as described in [3J: all reprocessing regions 
are centered on the coordinate origin and the wind axis is tilted inside the yz-plane leaning 
towards the positive y-axis. The observer's inclination and the tilting angle of the outflows are 
measured from the z-axis. The azimuthal viewing angle, </>, is taken inside the xy-plane with 
respect to the negative y-axis. The modeling is carried out using the latest version of the Monte 
Carlo radiative transfer code stokes presented in [5] and recently upgraded (Marin et al., in 
preparation) . 

We investigate the complex radiative coupling for a model composed of : (1) an optically 
and geometrically thick dusty torus, (2) an hourglass shaped, bi-conical polar wind filled with 
electrons, and (3) an equatorial, electron filled, radiation supported, flared disk. The equatorial 
scattering region is important in order to explain the observed polarization dichotomy between 
type-1 and type-2 AGN: it turns out from the observations that in type-1 AGN the polarization 
position angle in the optical/UV preferentially aligns with the projected symmetry axis of the 
polar outflows [6]. Observationally, this axis is assumed to be traced by the small-scale radio jets 
(note that in radio-quiet AGN these jets are much smaller and weaker than the ballistic, kpc- 
size jets seen in radio- loud objects). The small jets in radio-quiet AGN should indicate nuclear 
outflows that contribute to the larger-scale, bi-conical ionization regions seen in optical/UV 
imaging. When investigating the scattering-induced polarization produced by the different 
reprocessing regions listed above, it turns out that it is difficult to reproduce the observed 
aligned polarization at type-1 viewing- angles, unless an equatorial scattering region, such as a 
flared disk, is present (TJ El |9] . 



3. Spinning around NGC 1068 : modeling results 

The main goal in this research note is to revisit the X-ray polarimetry modeling of [3] while 
shifting the energy band to UV and optical wavelengths and without doing any morphological 
changes to the model (with the exception of leaving aside the lamp-post irradiation geometry 
used only in the X-ray modeling). 

3.1. Spectropolarimetric signature 

The spectropolarimetric modeling results are shown in Fig|TJ We consider four inclinations at 
which we investigate the spectral and polarization response of the system. The angle i ~ 26° 
features a type-1 object (face-on view), i ~ 60° refers to a viewing angle passing through the 
upper layers of the torus, i ~ 73° is around the expected inclination of NGC 1068 [TO] and, finally, 
i ~ 84° features an extreme type-2 object (edge-on view). We cover half a round in azimuthal 
angle at <fi = 35°, <fi = 80°, 4> = 125° and <f> = 170°. We extend the coverage to a full round for the 
spectra of the polarization position angle if). We define tp = 0° as "perpendicular" polarization 
with respect to the z-axis. It rotates clock-wise until ip = 90° ("parallel" polarization). 

Similarly to the X-ray results presented in [3] , the spectral flux is rather (^-independent at face- 
on view because the model does not deviate much from axis-symmetry. Since at low inclination 
the primary source is directly visible, the total flux is dominated by the power law continuum 
spectrum. Nevertheless, an additional scattered component contributes by more than 20% to 
the normalized flux. The flux level decreases rapidly as the observer's line of sight changes 
from type-1 to type-2 viewing directions. Heavy obscuration arises when the line of sight passes 
through the dusty torus, and the central source is no longer visible. Instead, scattered radiation 
dominates the spectrum. 

The spectral shape of the polarization percentage, P, depends on the inclination: at i ~ 26° 
and i ~ 60° , P has low values due to strong dilution by the unpolarized primary radiation coming 
from the directly visible source. The low polarization is nearly wavelength-independent in the 
UV, which reveals electron scattering inside the flared disk as the main polarization mechanism. 
The polar outflows are seen either in transmission (forward-scattering in the near cone) or in 
straight-line reflection (backscattering in the far cone) and they thus do not contribute much 
to the polarization. A dim contribution of dust scattering to the polarization appears toward 
increasing wavelengths when the Mie scattering phase function becomes less anisotropic so that 
multiply scattered photons have a higher probability to escape from the torus funnel. When 
the line of sight crosses the torus horizon and moves towards the highest inclination, the effects 
of dust scattering are less prominent as the torus funnel becomes more and more obscured. 
Electron scattering in the polar outflows is dominant and therefore P rises with i and remains 
constant with wavelength. This is in agreement with the observations [11, 12, 13J. The spectral 
shape of P is nearly independent of azimuth, but not its normalization, which varies with c/>. 

The polarization position angle ip shows a systematic behavior as the virtual observer travels 
around the object in the (^-direction. Despite the tilt of the outflows, the modeled ip confirms 
the observed polarization dichotomy at all ^-angles: for type-1 objects, we obtain ~ 90° and 
therefore the polarization roughly aligns with the axis of the outflow. For type-2 objects, the 
polarization is oriented perpendicularly to this axis. This picture would change, though, for 
more extreme tilting angles of the wind. The parallel polarization at type-1 viewing angles can 
only be preserved when the outflows are seen along a line of sight that mainly favors forward and 
backscattering. The polarization angle is symmetric in eft, which is expected because the axis of 
the inclined outflows lies in the yz-plane. The variation of ip with <ft is broader at type-2 than 
at type-1 view. By comparing the sign of ip at azimuthal viewing angles that are symmetric 
with respect to the yz-plane, we can recover the wind tilting direction: If the polarization angle 
is positive (negative), the polar outflows should be inclined in a clock- wise (counter clock- wise) 
direction from the z-axis. 
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Figure 1. Spectral and polarization results in the optical/UV for the NGC 1068 model with 
inclined polar outflows. From top to bottom the panels represent the four inclinations % ~ 26°, 
i ~ 60°, i ~ 73° and i ~ 84°. Left: spectral flux F/F cen t normalized to the source flux F cen t. 
Middle: percentage of polarization P. Right: polarization position angle tp. The different curves 
denote four azimuthal viewing angles (half a round) measured from the negative y-axis : <j) = 
35° (blue), <p = 80° (green), 4> = 125° (red) and cj) = 170° (black). For the polarization position 
angle, tp, the spectra at the symmetric azimuth (cj) = -cj)) are also shown using dashed lines. For 
ip = 0° the polarization is perpendicular to the z-axis; it rotates clock-wise with rising ip. 



The position of the outflow axis in NGC 1068 was investigated by [H] conducting kinematic 
modeling of the narrow line region. This region is assumed to be the natural continuation of 
the ionized winds. The kinematic modeling suggests that the axis does not lie in the plane of 
the sky but is indeed azimuthally rotated. From spectroscopy alone, however, this conclusion is 
difficult to verify and the polarization modeling carried out here can give an additional tracer 
of the scattering geometry. The combination of the polarization percentage and position angle 
adds two more observables that are sensitive to the position of the wind axis. 



3.2. Polarization imaging 

Results of our simulated polarization cartography are given in Fig. [2] We show the maps at an 
inclination of i ~ 73°, which is in agreement with the expected viewing angle towards NGC 1068 
as derived by [10J. The maps simultaneously show the polarized flux, PF V , and the polarization 
position angle, ip. The angle ip is represented by black bars drawn in the center of each spatial 
bin. A vertical bar indicates a polarization of ip = 90°, a bar leaning to the right denotes 
90° > ip > 0° and a horizontal bar stands for ■0 = 0°. For each pixel, the Stokes parameters are 
integrated over the full wavelength range. 

The images illustrate how the net polarization of the AGN is determined by the polarized 
flux integrated over the projected extension of the object. Independently of the azimuthal angle, 
strong polarized flux with a polarization position angle geared towards a vertical direction is 
spotted at the center. The polarization here is dominated by the combined scattering in the 
flared electron disk and the funnel of the relatively flat torus. It was shown in [5 J that both 
such scattering regions produce parallel polarization. Even though the central polarized flux 
predominates over the average polarized flux from the outflows by a factor of 5, the large spatial 
extension of the latter determines the net polarization. 

It is instructive to compare our model images with the observed UV polarization maps shown 
in Fig. 2 of [15]. Our model reproduces the arch-like distribution of polarization position angles 
across the extended outflows rather well. In the center, where according to |15j the main 
scattering source is located, the observed behavior still differs from our modeling. Although 
the observed polarization position angle near the central source shows some deviation from its 
behavior in the surroundings, this local turn is less prominent than our model suggests. It is 
possible that foreground matter on the observer's line of sight leads to additional polarization 
and can reduce the difference. Therefore, our analysis and modeling of the observed polarization 
images of NGC 1068 still needs some refinement, in particular in the central region of the object. 
It may also be that the inner surfaces of the torus are indeed steeper than assumed here so that 
the resulting polarization from the torus funnel is perpendicular and not parallel. Also, we have 
not investigated the impact of the electron optical depth in the polar outflows or the equatorial 
scattering disk yet. Such investigations are currently carried out and will be presented in a 
forthcoming paper (Marin et al., in preparation). 
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Figure 2. Modeling the spatial distribution of the polarized flux, PF U , of our model for 
NGC 1068 with inclined polar outflows. The PF V is denoted by the color-coded grid and 
normalized to the total flux. The polarization angle ip is represented by black bars (see text). 
The model region is observed at the inclination i ~ 73° for the same (^-angles as used in FigjTj 
Top, left: (j) = 35°, Top, right: (f> = 80°, Bottom, left: <j) = 125°, Bottom, right: 4> = 170°. 



